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Eli Waxman * ^^ ^ 

'^Dept. of Condensed Matter Physics, Weizmann Institute of Science, Rehovot 76100, Israel 



o 
o 
o 

(N 

Oh' 



> 

(N 

in 

a^ 
o 
o 
o 

Oh: 



X 



Observations suggest that 7-ray bursts (GRBs) are produced by the dissipation of the kinetic energy of a 
relativistic fireball. In this talk, recent work on the production of high energy neutrinos by GRB fireballs is 
reviewed. A significant fraction of GRB energy is expected to be converted to an accompanying burst of high 
energy neutrinos. Photomeson interactions produce a burst of ~ 100 TeV neutrinos in coincidence with the GRB, 
and a burst of ~ 10^* eV neutrinos following the GRB on a time scale of 10 s. Inelastic p-n nuclear collisions 
result in the production of a burst of ~ 10 GeV neutrinos in coincidence with the GRB. Planned 1 km^ neutrino 
telescopes are expected to detect tens of 100 TeV neutrino events, and several 10^* eV events, correlated with 
GRBs per year. A suitably densely spaced detector may allow the detection of several 10 GeV events per year. 
The detection of high-energy neutrino events correlated with GRBs will allow to constrain GRB progenitor models 
and to test the suggestion that GRBs accelerate protons to ~ 10^" eV. Moreover, such detection will allow to test 
for neutrino properties, e.g. flavor oscillations (for which upward moving r's would be a unique signature) and 
coupling to gravity, with an accuracy many orders of magnitude better than is currently possible. 



1. Introduction 

The widely accepted interpretation of the phe- 
nomenology of GRBs, bursts of 0.1 MeV-1 MeV 
photons lasting for a few seconds ||l|, is that the 
observable effects are due to the dissipation of 
the kinetic energy of a relativistically expanding 
wind, a "fireball," whose primal cause is not yet 
known (see pj for review). During the past three 
years, delayed low energy emission (afterglow) of 
GRBs has been discovered in X-ray, optical, and 
radio wave bands (see [0 for review). Afterglow 
observations confirmed the cosmological origin of 
the bursts, through the redshift determination of 
several GRB host-galaxies, and confirmed stan- 
dard model predictions of afterglows that result 
from the collision of an expanding fireball with its 
surrounding medium (see y for review). 

The physical conditions in the fireball dissipa- 
tion region allow Fermi shock acceleration of pro- 
tons to energy > 10^'^eV gP. The average rate 
at which energy is emitted as 7-rays by GRBs 
is comparable to the energy generation rate of 
> 10^^ eV cosmic-rays in a model where such 
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ultra-high energy cosmic-rays (UHECRs) are pro- 
duced by a cosmological distribution of sources 
1^,01 • These two facts suggest that GRBs and 
UHECRs have a common origin (see g for re- 
view). 

Emission of high energy neutrinos is expected 
to accompany the observed GRB photons, due to 
photomeson interactions of high energy protons 
with fireball photons y,|lO|_, and due to proton- 
neutron inelastic collisions |ll|-[l^. In this review, 
recent work on the production of high energy neu- 
trinos by GRB fireballs is summarized, and the 
implications for planned ^ 1km neutrino tele- 
scopes (the ICECUBE extension of AMANDA, 
ANTARES, NESTOR; see |l| for review) are dis- 
cussed. The fireball model is briefly described in 
§2. Proton acceleration in GRB fireballs is dis- 
cussed in §3. High energy neutrino production in 
fireballs and its implications for future high en- 
ergy neutrino detectors are discussed in §4. 

2. The fireball model 

In the fireball model of GRBs [|l5[, a compact 
source, of linear scale rg ^ 10^ cm, produces 
a wind characterized by an average luminosity 
L ~ lO^^ergs^^ and mass loss rate M = L/rjc^. 
At small radius, the wind bulk Lorcntz factor, T, 



grows linearly with radius, until most of the wind 
energy is converted to kinetic energy and T sat- 
urates at r ~ 7y ~ 300. Variability of the source 
on time scale At, resulting in fluctuations in the 
wind bulk Lorentz factor T on similar time scale, 
then leads to internal shocks [Q in the expand- 
ing fireball at a radius n « T^cAt. If the Lorentz 
factor variability within the wind is significant, 
internal shocks reconvert a substantial part of the 
kinetic energy to internal energy. It is assumed 
that this energy is then radiated as 7-rays by syn- 
chrotron and inverse-Conipton emission of shock- 
accelerated electrons. 

In this model, the observed 7-ray variability 
time, r^ Ti/V'^c ~ At, reflects the variability time 
of the underlying source, and the GRB duration, 
T ^ 10s, reflects the duration over which energy 
is emitted from the source. A large fraction of 
bursts detected by BATSE show variability on 
the shortest resolved time scale, '^ 10 ms p7| , 
and some show variability on shorter time scales, 
^ 1 ms [Q. This sets the constraint on under- 
lying source size, vq < cAt ^ 10^ cm. The wind 
must be expanding relativistically, with a Lorentz 
factor r ~ 300, in order that the fireball pair- 
production optical depth be small for observed 
high energy, - 100 MeV, GRB photons ||l|l. 

As the fireball expands, it drives a relativis- 
tic shock (blast-wave) into the surrounding gas. 
At early time, t <^ T, the fireball is little af- 
fected by this external interaction. At late time, 
t ^ T, most of the fireball energy is transferred 
to the surrounding gas, and the flow approaches 
the Blandford-McKee self-similar flow |§. The 
shock driven into the ambient medium continu- 
ously heats fresh gas, and accelerates relativis- 
tic electrons which produce through synchrotron 
emission the delayed radiation, "afterglow," ob- 
served on time scales of days to months. As the 
shock-wave decelerates, due to accumulation of 
ambient gas mass, the emission shifts with time 
to lower frequencies. For expansion into uniform 
density gas, the shock Lorentz factor decreases 
with time according to TBAiit) oc t^'^/*. 

During the transition to self-similar expansion, 
which takes place on (observed) time scale com- 
parable to T, reverse shocks propagate into the 
fireball ejecta and decelerate it. At this stage 



the shocked plasma expands with the self-similar 
Lorentz factor, which for expansion into uniform 
density gas is given by (e.g. |10|) 
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while the unshocked fireball ejecta propagate with 
the original expansion Lorentz factor, F = 77 > 
TBAiit = T). Here, E = lO'^'^E^:^ erg is the total 
fireball energy, T — lOTi s, and n = \ riocm^"^ is 
the ambient gas density, no = 1 is typical to the 
interstellar medium. 

Internal shocks are expected to be "mildly" rel- 
ativistic in the fireball comoving frame, i.e. char- 
acterized by Lorentz factor F^ — 1 ~ 1 in the 
wind rest frame, since adjacent shells within the 
wind are expected to expand with Lorentz fac- 
tors which do not differ by more than an order of 
magnitude. Moreover, the reverse shocks are also 
expected to be mildly relativistic, since the ratio 
v/^BMit = T) is not far from unity. 

3. UHECRs from GRB fireballs 

3.1. Fermi acceleration in GRBs 

In the fireball model, the observed GRB and 
afterglow radiation is produced by synchrotron 
emission of shock accelerated electrons. In the 
region where electrons are accelerated, protons 
are also expected to be shock accelerated. This 
is similar to what is thought to occur in super- 
novae remnant shocks [|l[. We consider below 
proton acceleration in internal and reverse fire- 
ball shocks. Since these shocks are mildly rela- 
tivistic, we expect results related to particle ac- 
celeration in sub-relativistic shocks (see, e.g., ||2l|] 
for review) to be valid for the present scenario. 
In particular, the predicted energy distribution 
of accelerated protons is dNp/dcp oc e~^. 

Two constraints must be satisfied by fireball 
wind parameters in order to allow proton accel- 
eration to e„ > 10^" eV in internal shocks [Hi: 



Cs/6 > 0.02F|5e2 2^L-i2, 
and 



F > imel%M-T- 



(2) 



(3) 



Here, e^ 



102°ep,2o eV, At = lO-^At.z s, T 



10^'^r2.5 is the plasma expansion Lorentz factor, 
Lj — 10^^L-y_52erg/s is the 7-ray luminosity, £^b 
is the fraction of the wind energy density which 
is carried by magnetic field, 47rr^cr^(B^/87r) — 
^bL, and ^e is the fraction of wind energy car- 
ried by shock accelerated electrons. To obtain 
the constraints for acceleration in reverse shocks, 
Ai should be replaced by T. Since the electron 
energy is lost radiatively, L^ k, ^^L. The first con- 
straint must be satisfied in order for the proton 
acceleration time ta to be smaller than the wind 
expansion time. The second constraint must be 
satisfied in order for the synchrotron energy loss 
time of the proton to be larger than ta- 

The constraints that must be satisfied to allow 
acceleration of protons to energy > 10^" eV are 
remarkably similar to those inferred from 7-ray 
observations. F > 100 is implied by observed 7- 
ray spectra, and magnetic field close to equipar- 
tition, ^B ~ 1, is required in order for electron 
synchrotron emission to account for the observed 
radiation. 

We have assumed in the discussion so far that 
the fireball is spherically symmetric. However, 
since a jet-like fireball behaves as if it were a 
conical section of a spherical fireball as long as 
the jet opening angle is larger than F"^, our re- 
sults apply also for a jet-like fireball (we are inter- 
ested only in processes that occur when the wind 
is ultra-relativistic, F ~ 300, prior to significant 
fireball deceleration). For a jet-like wind, L in 
our equations should be understood as the lumi- 
nosity the fireball would have carried had it been 
spherically symmetric. 

3.2. Energy generation rate 

The observed GRB rcdshift distribution im- 
plies a typical GRB 7-ray energy oi E r^ 10^'^ erg, 
and a GRB rate of i?GRB ~ S/Gpc'^yr at z ^ 1. 
The present, z = 0, rate is less well constrained, 
since most observed GRBs originate at redshifts 
1 < 2; < 2.5 Q. Present data are consistent with 
both no evolution of GRB rate with redshift, and 
with strong evolution (following, e.g., star forma- 
tion rate), in which i?GRB(z — 1)/Rgkb{z = 0) '^ 
10 |24|. The energy observed in 7-rays refiect the 
fireball energy in accelerated electrons. If accel- 
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Figure 1. The UHECR flux expected in a 

cosmological model, where high-energy protons 
are produced at a rate {eidhp / dep) z=a = 0.8 x 
lO^^'erg/Mpc yr as predicted in the GRB model 
[Eq. (^)], compared to the Fly's Eye, Yakutsk 
and AGASA data, la flux error bars are shown. 
The highest energy points are derived assuming 
the detected events (1 for Fly's Eye and Yakutsk, 
4 for AGASA) represent a uniform flux over the 
energy range 10^° eV-3 x 10^" eV. 



erated electrons and protons carry similar energy 
(as indicated by afterglow observations |2^) then 
the GRB cosmic-ray production rate is 



elidnp/dep),^o = lO^^Cerg/MpcVr, 



(4) 



with C in the range of ~ 10-°-5 to - 10°-5. This 
result is consistent with more detailed analyses 
pq ], where RQusiz = 0) x £■ in the range of 0.2 
to 2 times 10^'*erg/Mpc yr is obtained (under the 
assumption that -Rgrb is proportional to star for- 
mation rate) for different models of star formation 
rate evolution. 

In Fig. 1 we compare the UHECR spectrum, 
reported by the Fly's Eye I?), the Yakutsk g, 
and the AGASA ||2^ experiments, with that pre- 
dicted by the GRB model (assuming GRB rate 
follows star formation rate). The generation rate 
(0) of high energy protons is remarkably simi- 



lar to that required to account for the flux of 
> 10^® eV cosmic-rays (The flux at lower ener- 
gies is most flkely dominated by heavy nuclei of 
Galactic origin [2^] , as indicated by the flattening 
of the spectrum at « 10^^ eV and by the evidence 
for a change in composition at this energy |3C|]). 

The suppression of model flux above 10^^'^ eV 
is due to energy loss of high energy protons in in- 
teraction with the microwave background, i.e. to 
the "GZK cutoff" [|l|. The available data do not 
allow to determine the existence (or absence) of 
the "cutoff" with high confidence. The AGASA 
results show an excess (at a ^ 2.5a confidence 
level) of events compared to model predictions 
above 10^'^eV. This excess is not confirmed, how- 
ever, by the other experiments. Moreover, since 
the lO^^eV flux is dominated by sources at dis- 
tances < 100 Mpc, over which the distribution of 
known astrophysical systems (e.g. galaxies, clus- 
ters of galaxies) is inhomogeneous, significant de- 
viations from model predictions presented in Fig. 
1 for a uniform source distribution are expected at 
this energy |0] . Clustering of cosmic-ray sources 
leads to a standard deviation, a, in the expected 
number, N, of events above 10^° eV, given by 
a/N = 0.9(do/10Mpc)°-^ ||, where do is the un- 
known scale length of the source correlation func- 
tion and do ~ 10 Mpc for field galaxies. 

Although the rate of GRBs out to a distance of 
100 Mpc from Earth, the maximum distance trav- 
eled by > lO^'' eV protons, is in the range of 10^^ 
to 10~'^yr^^, the number of different GRBs con- 
tributing to the fiux of > 10^° eV protons at any 
given time may be large. This is due to the dis- 
persion T in proton arrival time, which is expected 
due to deflection by inter-galactic magnetic fields 
and may be as large as 10^ yr, implying that the 
number of sources contributing to the flux at any 
given time may be as large as t x lO^^yr^^ = 10'^ 



4. High energy Neutrinos 

4.1. Photomeson production 
4.1.1. Internal shocks neutrinos, ^ 10^^ eV 
The decay of charged pions, 7r+ — > /i+ + Vfj, —^ 
e+ + i^e + I7p + j/p, created in interactions be- 
tween fireball photons and accelerated protons. 



results in the production of high energy neutri- 
nos [g[. The key relation is between the observed 
photon energy, e-y, and the accelerated proton's 
energy, e^, at the photomeson threshold of the 
A-resonance. In the observer frame. 



e^ep = G.2GcV^r^ 



For F 



300 and e^ 



(5) 

1 MeV, the typical 
observed energy of photons emitted in internal 
shocks, we see that characteristic proton ener- 
gies ^ 10^^ eV are required to produce neutrinos 
from pion decay. Typically, the neutrinos receive 
^ 5% of the proton energy, leading to neutrinos 
of - IQi"* eV. 

The fraction /^(e^) of proton energy lost to 
pion production is given by the ratio of the wind 
expansion time, FAi in the wind frame at the in- 
ternal shock stage, and the proton photopion en- 
ergy loss time. For observed GRB photon spec- 
tra, dNj/dej oc e~^ with /3 ~ 1 for photon en- 
ergies below the break energy e^ « 1 MeV and 
/? ~ 2 for e^ > e'^, f^ is independent of energy for 
Cp > 0.2F2GeVV4 ~ 10^^ eV, 



U « 0.2- 



^7,52 



eVlMeV)F|.5Ai_2- 



(6) 



and /tt oc Ep for Ep < 10^® eV. For a typical burst, 
E « 10^'^ erg at redshift z ^ 1, Eq. (||) implies a 
neutrino fluencc of |9| 



e^^.^. 
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where a = for e^ > lO^^'eV and a = 1 for 
El, < lO^'^eV, and Vx stands for i'^, Dp and i^e- 

The neutrino spectrum (|^) is modified at high 
energy, where neutrinos are produced by the de- 
cay of muons and pions whose life time exceeds 
the characteristic time for energy loss due to adia- 
batic expansion and synchrotron emission [pSlJsfl . 
The synchrotron loss time is determined by the 
energy density of the magnetic field in the wind 
rest frame. For the characteristic parameters of 
a GRB wind, synchrotron losses are the domi- 
nant effect, leading to strong suppression of u flux 
above - lO^*' eV. 



4.1.2. Afterglow neutrinos, ~ 10^^ eV 

During tlie transition to self-similarity, protons 
and electrons are accelerated to high energy in 
the reverse shocks. High energy protons may in- 
teract with the 10 eV-1 keV photons radiated by 
the electrons, to produce through pion decay a 
burst of duration '^ T of ultra-high energy, 10^^- 
10^^ eV, neutrinos as indicated by Eq. (||) [|0|. 

While afterglows have been detected in sev- 
eral cases, reverse shock emission has only been 
identified for GRB 990123 ||. Both the detec- 
tions and the non-detections are consistent with 
shocks occurring with typical model parameters 
|23| , suggesting that reverse shock emission may 
be common. The predicted neutrino emission 
depends, however, upon parameters of the sur- 
rounding medium that can only be estimated 
once more observations of the prompt optical af- 
terglow emission are available. 

If the density of gas surrounding the fireball is 
n ~ Icm"'^, a value typical to the inter-stellar 
medium and consistent with GRB 990123 obser- 
vations, then the synchrotron emission of reverse 
shock electrons is expected to peak in the X- 
ray band, e!j, « 1 keV, with luminosity Lx ~ 
2 X 10^°erg/s |1^. Using these parameters in Eq. 
(^, replacing At with T « 10 s and recalling 
that r = 250 at the reverse shock stage (see §2), 
we find /tt « 0.01 for protons of energy 10^^ eV. 
Thus, the expected neutrino fluence for a typical 
burst, E « 10^3 grg at z ~ 1, is [H 
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(8) 



where a = 1/2 for e^ > lO^^eV and a = 1 for 
e^ < lO^^eV. Here too, v^ stands for f^, v^ and 
v^. The value of a, a = 1/2 for e^ > 10^'^ eV, 
corresponding to photomeson interactions of pho- 
tons with energy below the break e!j, « 1 keV, dif- 
fers from the value a = in the case of internal 
shocks discussed in §4.1.1, since the low energy, 
€j < e^, photon spectrum is different in the two 
cases: dN^/de-y oc e~^ with /? = 1 for internal 
shocks, and (3 — 3/2 for reverse shocks. 

Some GRBs may result from the collapse of a 
massive star (e.g. [Q), in which case the fireball 
is expected to expand into a pre-existing wind. 
For typical wind parameters, the transition to 



self-similar behavior takes place at a radius where 
the wind density is n m lO^cm^^ ^ lcm~^. 
The higher density implies a lower Lorenz fac- 
tor of the expanding plasma during the transi- 
tion stage, and hence a larger fraction of proton 
energy lost to pion production. Protons of energy 
Cp > 10^* eV lose all their energy to pion produc- 
tion in this case, and a typical GRB at z '-~^ 1 is 
expected to produce a neutrino fluence |lQ,pq] 






" GeV 



cm^ 



(9) 



where a = for e^y > lO^^eV and a = 1 for 
e^ < IQi^eV. 

The neutrino flux is expected to be strongly 
suppressed at energy > 10"^^ eV, since protons 
are not expected to be accelerated to energy ^ 
1020 eV. 

4.2. Inelastic p-n collisions 

The acceleration, F oc r, of fireball plasma 
emitted from the source of radius tq (see §2) is 
driven by radiation pressure. Fireball protons 
are accelerated through their coupling to the elec- 
trons, which are coupled to fireball photons. Fire- 
ball neutrons, which are expected to exist in most 
progenitor scenarios, are coupled to protons by 
nuclear scattering as long as the comoving p- 
n scattering time is shorter than the comoving 
wind expansion time r/Tc = ro/c. As the fireball 
plasma expands and accelerates, the proton den- 
sity decreases, Up oc r~^T~^, and neutrons may 
become decoupled. For rj > rjpn, where 

r]pn « 400Ll!^'^rQ^^^ 



(10) 



and To = 10 ro.7 cm, neutrons decouple from the 
accelerating plasma prior to saturation, T — rj, 
at F = rjpii ?7~^/^ [Q. In this case, relativistic 
relative velocities between protons and neutrons 
arise, which lead to pion production through in- 
elastic nuclear collisions. Since decoupling occurs 
at a radius where the collision time is similar to 
wind expansion time, each n leads on average to 
one pair of iyf>. The typical comoving neutrino 
energy, ^ 50 MeV, implies an observed energy 
^ 10 GeV. A typical burst, E = 10^^ erg at 
z ~ 1, with significant neutron to proton ra- 
tio and 1] > 400 will therefore produce a flu- 



ence F{i^e + v^) ~ {).hF{v^ + v^) ^ lO^^'cm^^ 
of ~ 10 GeV neutrinos. 

Relativistic relative p-n velocities, leading to 
neutrino production through inelastic colhsions, 
may also result from diffusion of neutrons be- 
tween regions of the firebaU wind with large dif- 
ference in r [y_3|. If, for example, plasma expand- 
ing with very high Lorentz factor, F > 100, is 
confined to a narrow jet surrounded by a slower, 
r ~ 10 wind, internal coUisions within the slower 
wind can heat neutrons to relativistic tempera- 
ture, leading to significant diffusion of neutrons 
from the slower wind into the faster jet. Such 
process may operate for winds with 77 < 400 as 
well as for rj > 400, and may lead, for certain 
(reasonable) wind parameter values, to ~ 10 GeV 
neutrino flux similar to that due to p-n decoupling 
in a ?7 > 400 wind. 

4.3. Implications 

The predicted fluence of ~ 10^^ eV neutrinos 
produced by photomeson interactions in internal 
fireball shocks, Eq. (0), implies a detection prob- 
ability ~ 10~^'^ per burst in planned 1km neu- 
trino telescopes |1J], corresponding to detection 
of several tens of muon induced neutrinos per 
year correlated in time and direction with GRBs, 
given the observed GRB rate of « 10"^yr~^. The 
predicted fluence of ~ 10^^ eV neutrinos, pro- 
duced by photomeson interactions during the on- 
set of fireball interaction with its surrounding 
medium, Eqs. (||,||), implies a detection probabil- 
ity r^ lO""''^ in the case of fireball expansion into 
typical inter-stellar medium gas, and '^ 10~^'^ in 
the case of fireball expansion into a pre-existing 
massive star wind. Several muon induced neutri- 
nos per year in a 1km detector are expected in 
the latter case. 

Detection of high energy neutrinos will test the 
shock acceleration mechanism and the suggestion 
that GRBs are the sources of ultra-high energy 
protons, since > 10^^ eV (> 10^* eV) neutrino 
production requires protons of energy > 10^^ eV 
(> 10^9 eV). The dependence of - lO^-* eV neu- 
trino flux on wind Lorentz factor, Eq. (|6|), and 
the dependence of ~ 10^^ eV neutrino flux on fire- 
ball environment, imply that the detection of high 
energy neutrinos will also provide constraints on 



wind Lorentz factor ||3^, and on the GRB pro- 
genitor. 

Inelastic p-n collisions may produce ~ 10 GeV 
neutrinos with a fiuence of ~ 10~^cm~^ per burst, 
due to either p-n decoupling in a wind with high 
neutron fraction and high, > 400, Lorentz factor 
||ll| , |l2|| , or to neutron diffusion in a wind with, 
e.g., strong deviation from spherical symmetry 
p^ . The predicted number of events in a Ikm^ 
neutrino telescope is ~ 10yr~^. Such events may 
be detectable in a suitably densely spaced detec- 
tor. Detection of ^ 10 GeV neutrinos will con- 
strain the fireball neutron fraction, and hence the 
GRB progenitor. 

Detection of neutrinos from GRBs could be 
used to test the simultaneity of neutrino and pho- 
ton arrival to an accuracy of ^ Is (^ 1 ms for 
short bursts) , checking the assumption of special 
relativity that photons and neutrinos have the 
same limiting speed. These observations would 
also test the weak equivalence principle, accord- 
ing to which photons and neutrinos should suffer 
the same time delay as they pass through a grav- 
itational potential. With 1 s accuracy, a burst at 
100 Mpc would reveal a fractional difference in 
limiting speed of 10~^^, and a fractional difference 
in gravitational time delay of order 10~^ (con- 
sidering the Galactic potential alone). Previous 
applications of these ideas to supernova 1987A 
(see Ig^ for review), where simultaneity could 
be checked only to an accuracy of order several 
hours, yielded much weaker upper limits: of or- 
der 10~® and 10~^ for fractional differences in the 
limiting speed and time delay respectively. 

The model discussed above predicts the pro- 
duction of high energy muon and electron neu- 
trinos. However, if the atmospheric neutrino 
anomaly has the explanation it is usually given, 
oscillation to i^t's with mass ~ 0.1 eV BS], then 
one should detect equal numbers of v^'s and Vr^s. 
Up-going r's, rather than /i's, would be a dis- 
tinctive signature of such oscillations. Since i^r^s 
are not expected to be produced in the fireball, 
looking for t's would be an "appearance experi- 
ment." To allow flavor change, the difference in 
squared neutrino masses, A?n^, should exceed a 
minimum value proportional to the ratio of source 
distance and neutrino energy |3^. A burst at 



100 Mpc producing 10^'^eV neutrinos can test for 
Am? > 10~^^eV^, 5 orders of magnitude more 
sensitive than solar neutrinos. 
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